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Abstract 
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the method of inhomogeneity expansion. We first treat the non-relativistic case 
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we also add a 4-fermi type interaction and show that under certain circumstances, 
it corresponds to a Zeeman type term in the effective action. 
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1. Introduction 

The behavior of a system of interacting planar fermions in an external magnetic field 
has many interesting features. In 2+1 dimensions such systems are of relevance for the 
description of the Quantum Hall Effect and high T c superconductivity. In addition, we 
have the ulterior motive of adapting our technique to 3+1 dimensional situations. It is 
of particular (and potentially of practical) interest to study Abelian (or non-Abelian) 
theories involving fermions (in 3+1 dimensions) in the environment of strong magnetic 
fields. This is especially rewarding in a cosmological context, where strong magnetic fields 
may be present in the early universe, collapsing supernovae and perhaps superconducting 
strings. To give an idea of the magnitudes of the magnetic fields involved, we may 
say that those associated with compact astrophysical objects range between O(10 4 ) T 
(magnetic white dwarfs) and O(10 10 ) T (supernovae) [|l]]. In between, neutron stars have 
associated magnetic fields of O(10 8 ) T. It has also been suggested in some models for 
extragalactic gamma bursts that very strong magnetic fields O(10 13 ) T are involved 0. 
The vacuum structure of gauge theories in a strong magnetic background, has been and 
continues to be a field of considerable interest 0,0. The minimal standard model also 
exhibits interesting features in the presence of a strong magnetic background, through the 
self-energy corrections of the neutrinos, through vacuum polarization and also through 
non-perturbative fermion number violation ||, [^J. It has also been observed that 
in certain planar fermionic theories with four-fermi interactions, a strong magnetic field, 
in conjunction with a finite density, induces breaking of the chiral symmetry and thereby 
generating dynamical fermion mass. We can address all these issues through minimal 
modifications of the techniques described below. The effective action in a magnetic 
background plays a pivotal role in all of them. This is one of our primary motivations 
in outlining our method in considerable detail. 

Effective actions neatly encapsulate the contributions of quantum corrections to the 
classical action and thus permit a classical analysis of fully quantum problems, an advan- 
tage not to be lightly regarded. Given an effective action for a quantum theory, one can 
read off the various physical quantities of interest like effective masses, permeability and 
susceptibility with relative ease. In this paper, we wish to find the effective action of the 
system in terms of the fluctuating gauge fields that the planar fermions are subjected to, 
by tracing the fermions out of the partition function. The formal process of integrating 
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Grassman valued fields out yields a functional determinant, generically a nonlocal quan- 
tity. In order to obtain a Lagrangian comprising local fields and their derivatives, we 
have to expand the determinant within some controlled scheme. The derivative expan- 
sion method, also called the inhomogeneity expansion in the present context, is one such 
scheme. To implement this, we assume that the fluctuating fields coupled to the fermions 
vary slowly in space-time compared to the natural scales set by the background magnetic 
field. Namely, the typical momenta associated with these fluctuations are <C j, where 
/ = -7= is the magnetic length and the typical frequency is to <C uj c , where u) c m is the 
cyclotron frequency. This naturally sets a limit on contributions from the higher deriva- 
tives of the electromagnetic fields by filtering out and discarding the high frequency and 
the high momentum modes and thus enables us to write a local Lagrangian containing 
the desired number of derivatives of the field. 

Our method obviously works well in the absence of gaps in the spectrum and fails 
precisely at the points where such gaps exist. Particularly, in this problem, where the 
spectrum is described by Landau levels, the method expectedly fails where the gap 
between successive levels is probed. In practical terms, we see, in our subsequent cal- 
culations, the emergence of spurious singularities in the form of derivatives of the delta 
function precisely at these points. In this work, we have adopted the somewhat prag- 
matic stance of disregarding these singularities and tacitly admitting the shortcomings 
of our method at these particular points. We however hasten to add that away from 
these particular points, derivative expansion is a perfectly reliable method and leads to 
a rather formulation of the effective action. If thermal effects are taken into account 
resulting in the smearing out of these discrete levels through thermal fluctuations, this 
method can be readily applied to even the points mentioned earlier. 

Although the solution of the "free" single particle problem is known in this case, 
the construction of an effective action (or potential) for the interacting theory is rather 
involved. In particular since the theory is gauge invariant to begin with, one would like 
to find a method or rather an expansion, that will exhibit this gauge invariance order by 
order in a manifest way. 

In this paper we construct such a method. We apply it to the nonrelativistic test caseS 

J This problem was originally considered by Ray and Sakita Q for the non zero temperature case. 
Our results agree with theirs for the T=0 case (see also M) 
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before tackling the fully relativistic theory. The basic idea of our method is the use of 
guiding center coordinates which are the natural variables for this problem. Using these 
variables, we are able to convert the problem of calculating the fermion determinant into 
the problem of calculating matrix elements of known commutators and of lowering and 
raising operators. In addition we make several unitary transformations along the way to 
expose the inherent gauge invariance of the system at each stage of the approximation. 

The organization of this paper is as follows: we start with the nonrelativistic problem, 
establishing the notation and describing the approximations entailed and transformations 
made in section 2. In section 3 we present the calculations and the results for the nonrel- 
ativistic case. In section 4 we study the effect of adding a four Fermi type interaction to 
the original Lagrangian and show that in a nonrelativistic theory this kind of term, under 
certain circumstances, can induce a Zeeman type interaction in the effective Lagrangian. 
In section 5 we consider the fully relativistic QED . We conclude with discussions in 
section 6. Useful relations are deferred to the appendices. 

2. Nonrelativistic Q.E.D 

To set the stage for the calculation we start by establishing the notation. Consider the 
single particle Landau problem. The dynamics of a single planar particle in an uniform 
magnetic field normal to the plane is described by the Landau hamiltonian: 

H = — P 2 

° 2m 

where A may be chosen in a specific gauge A = B(—y,x). The wave functions of the 
system, (x\n,X), are labeled by two parameters; n and X. n, which is an integer is 
the Landau level index and X, a continuous parameter, measures the degeneracy of any 
given Landau level. The single particle energies are given by: 

E n = (n+-)w c (1) 

where w c = ^ is the cyclotron frequency. These energies are obviously independent of 
X and the operator X which measures the value of X commutes with the hamiltonian. 
We call X and its conjugate Y the guiding center coordinates. They play a central role 
in our subsequent discussions. 



4 



Let us define tli = pi + |x 2 and n 2 = ]3 2 — yXi where [n 1 ,n 2 ] = p. We can 
similarly define a = -^(Jli + «n 2 ). The guiding center operators are then defined as 

X = X! + / 2 fl 2 and f = x 2 -Z 2 fli (2) 

thus [X, Y\ = —il 2 . These are the two sets of cannonically conjugate operators that we 
need: 

a^a\n, X) = n\n, X) X\n, X) — X\n, X) 

The normalization for X is fixed by the degeneracy of the Landau levels per unit area. 
Namely 

JdX\(n,X\n,X)\ 2 = ^ f2 (3) 

-, -y- -y- 

where (x\n,X) = ^ e'i 7 "■» ra ( £ =— ) and u n (^j— ) are the standard harmonic oscilator 
wave functions. 

Using an euclidean metric, the partition function for the corresponding many particle 
system is given by: 

Z = J D^D^e~J d * dt ^^ + ^W 2 -^ (4) 

where Il r = — id T — iA = p T — eA T ,and II = p r — eA. Integrating the fermions out we 

get Z = e~ Se ff where 

S eff = -Trlog[in T + ^-tl 2 - fi] (5) 



The corresponding current is 



5S ef.f 



where J„(x) = (ie— t-^- — , {it,-, — t—^t- — !)• We translate the bra and ket to 

x — 0, consequently the operators x^ in J M get translated to x^ + x M and we expand all 
the functions of x^ + x^ around x M (derivative expansion). 

To make the gauge invariance as manifest as possible perform a succession of unitary 
transformations T = WVU, on J M where 

JJ(x^t} — e -ie(^-A+^it i it j did j A j + ^i<i i ^ji<i k a i d j d k A+...) ^ 

V(x,t) = e - ie irAr(x)+^-A T (x)+^-A T + ... 

W{x,t) = e^™ 
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We notice that T, acting on |x M = 0), produces just a factor of unityo. Thus 

(j,(x,t)) = (o|rJ,r t |o). (7) 

Since A T is entirely fluctuating, we refer to it as a T , similarly the fluctuating part of Aj 
is referee! to as a^. 

After all these transformations we obtain: 

(J T ) = *e(0\-„ 1 - s 10) (8) 

U) = —(oi in,, — R \\o) 

W 2m VI | l ' jffo + _L(n- A + A-II + A 2 )+*A T ' 
where if = «n T + - // , 

i?o|n, = r n |n,X,w) = [uy — //+ (n + -)w c ]|n,X, w) 

and 

A T = — e[x ■ 9a T + -x,;Xj9j9ja r + Xjfaja T + . . .] (9) 

A, = -e[fai + \r 2 'di + fx • <9d; + ^x^ + ^XjX,-^^ + . . .]. 
2 2 3 

i^o can be solved exactly (its spectrum being the L.L) and the rest is treated perturba- 
tively in the following section. In fact, 

H \n, X, w) = T n \n, X, w) = [iw — fi + (n + -)w c }\n, X, w) . (10) 
3. Effective Action 

The most efficient way of computing the effective action is to obtain (11^) the po- 
larization tensor, from (J M ) and write: 

S e ff = - J dxdya^x^IL^x - y))a„(y) (11) 

where (IT^) = ^H^) • Using the remaining rotational symmetry, the only independent 
components of (n^) that are needed to reconstruct the effective action are (n o), (n i), 



2 At this point we note that the transformation V fails to be unitary in the corresponding finite 
temperature problem due to the compactness of the domain of definition of r , (0 — > (3 = 1/T). At 
T = 0, however, V is unitary. 
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(rin) and (II12). In section 2 we transformed (J M ) to a form that allows the use of 
perturbation theory. Thus for example we can expand eq.(||) to get 

W ie{0\- 1 10) (12) 

ip T -/i + iA T + ^(n + A) 2 

= «<0|^|0> - ze(0|i- (zA r + ±-(fi ■ A + A ■ ft + A 2 )) ±-\0) + . . . 

Now the A^ generically contain We use the definitions of X and Y from section 2 
to convert the 2^ into X and F and Ilj. The action of X and Y on |x M ) = (or its 
complex conjugate) enables us to get rid of X and Y altogether and obtain an operator 
that depends only on II. Subsequently we can introduce a complete set {\n, X)} of 
states and integrate over X (as the operator does not contain any further guiding center 
coordinates) immediately to obtain the degeneracy factor p — 2^72- The remaining 
evaluation of harmonic oscillator matrix elements is straightforward. 

To clarify the foregoing statements through an example, we compute the contribution 
proportional to fu to (J T ). This is given by 

l^oii-^-xA^lo) 

now, X4 = X — / 2 n 2 and x 2 = Y + l 2 fli , therefore 

(o\x = / 2 (o|ri 2 (o\y = -z 2 (o|ri 2 . 

Also [n 1; i] = -^n 2 ^ and similarly for U 2 , [tl 2 ,^] = ^Hi^. Thus the 
above matrix element becomes: 

!^ ( o|[n 2 ,i-]n 2 i- + [n I ,i-]n 1 i-|o) (13) 
g2 (o| 1 ri 1 ri 1 1 n 1 n 1 |o) 

2m 2 H H H q H H H 



■^(oi^n^n^io) 

™ -"0 -no -"0 



4 £ /«/^ < '"^-;><";' X ' ro ' > K|n 1 M ( n 2 |n 2 |n3) 



ni,ri2,n3 " " 7"i]_ 722 723 

The integration over X can now be done trivially yielding ^p^m,^, reducing the matrix 
element to: 

e 2 f dw (ni| III I n 2 )(n 2 1 II2 1 n 3 ) 



^ / 27T r 2 r 
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Using the definition of II we get 

ie 2 r dw 00 



Anm 2 l A 

and we use 



r dw . 

/^§ n + 1) 



(14) 



j i_ 



n=0 

Thus the above matrix element reduces to (retaining 



r n r n +i 

only the non-singular part as discussed in the introduction) 

ip 2 00 f dw 1 ie 2 00 1 

— T/ ; tt — fi2(x, r) = — 6(u> c (n + -) — a) /Wx, r) 

27r^ i 2 7 rra;- /U +(n+|)u>/ m ; 2tt ^ 1 cl 2 ; WJm ' ; 

For operators involving f we use [f , -^-] = and proceed as above. Another simplifying 
feature is that any matrix element with an odd number of Ilj is zero. 

It is easy to show that all terms which could possibly lead to gauge noninvariance 
vanish identically. We append an example here. The term proportional to Oj in ( J T ) leads 
to a gauge noninvariant term in the effective action. The matrix element proportional 
to this term is 

^ (0 |_Ln^i-|o> ai + (i^2) (is) 

2 1 m\ _ 2 



Now since t 2 -^|0) = ^|0) we get 



ie 2 1 - 1 

— (0| — (IIA + I^) — 10> = (16) 
m H H Q 

so that in fact this term does not appear. 

At this point we collect all the relevant terms, i.e. those proportional to dfa T , d\d T ai 
and /12 for (J T ) and proportional to d 2 a±, d 2 a 2 , did T a T , d 2 d T a T and d 2 f± 2 for (Ji) and 
obtain : 

(Tlrr(x,y)) = -^C(d 2 l+ d 2 2 )5(x-y) + ... (17) 

Z7TW C 

p 2 ie 2 
(U Tl (x,y)) = -- Cd T d 1 S(x-y)-—Cd 2 S(x-y) + ... 

ZnW c LlX 

2 • 2 

P IP 

(U T2 (x,y)) = -- Cd T d 2 + — Cd 1 5(x-y) + ... 

2nw c In 
2 2 

(n n (x,y)) = J—C%8{x-v)-—Cx%8(x-y) + ... 
ZttWc Tim 

p 2 ip 2 
(U 12 (x,y)) = —C 1 d 1 d 2 5(x-y)-—Cd T 5(x-y) + ... 
Tim Zrr 

(U 22 (x,y)) = J!-C%6(x - y) - —CtfSix - y) + . . . 
2ttw c Txm 



where 

» f dw / 9 \ 1 ~ ^, , L N 



J> + 2» / 2? I 1 + re J ?: = £<» + 2 )e( " - (n + 2 >«■> + • • • • (19) 



and Ci is 
Ci 

The effective action is reconstructed as : 

S eff = | rf 3 xa M ^a p + ^-C j d 3 x £ 2 + -^—C x J d?x B 2 + . . . (20) 

where S and B are the euclidean fluctuating electric and magnetic field strength respec- 
tively. If we choose /i such that exactly N Landau Levels are filled C —>■ N + 1 and 

°1 2 • 

We note two points regarding this effective action. The first is the emergence of 
a C-S term. The second is that since relativistic invariance is broken by the external 
magnetic field, we get different coefficients for the fluctuating £ 2 and B 2 terms. Further 
if the magnetic field is very strong (i.e. I <C — ), the second term in the effective action 
does not contribute. This corresponds to dimensional reduction. We also observe that if 
there are no fermions ,i.e. // = 0, the magnetic B 2 term in the gauge field action is not 
renormalized by quantum effects. This is to be contrasted with the relativistic case (see 
eq.flSBD) where the vacuum (filled Dirac sea) does contribute to the renormalization of 
this term. 

We can use this effective action (rotated to minkowski space) to obtain physically 
interesting quantities, for example the dielectric function 



e(p) = l+ =1 + o C 



n 00 (o,p) , , e 2 

= i-t- 

and the magnetic permeability 

i _ t + nKo^p) = 1 _ ± Ci _ 



/i(0, p) p 2 rrm 

4. Effective Zeeman Interaction 

Having obtained the effective action we take this oportunity to explore an interesting 
phenomenon namely the effect of adding a four fermi term. In this section we show 
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that if a four fermi term is added to the original Lagrangian introduced in the previous 
sections, an effective Zeeman-type interaction could be induced in certain cases. 



The partition function for the system is 
Z = J D^DipD^e^'^ 



(21) 



% \ cTx 



^[ido -eA + fx + — (iV + eA) 2 ]^ - g(^^) 2 
Using an auxiliary field <fi to disentangle the quartic term, we obtain 



Z = J Di/jD^exp ji J S 



x 



^[id - eA + n + ^-(iV + eA) 2 + fyajpj) + i 



(22) 

We wish to show that the effective potential for obtained upon integrating the fermions 
out is minimized by a non-zero constant configuration which we call 0o- This in turn 
would mean that the Yukawa term in the action above yields an effective Zeeman term 

The effective potential for <j> is given by 



Veff = 



1 



2tt/ 2 ^ 7 2m 



dki) 



A; - £■„ + /x + 



6 1 

2 



(23) 



where _E n = ^j?(n + |) and b = \^2g<f). The presence of the sum of two logarithms is due 
to the two component nature of the fermion. The mean density is computed from V e ff 
as (0|p|0) = -f^, which yields 



+ 



1 



k - £ n + /i + | k - E n + [i- 



(24) 



On performing the &o integral, we get 



(0|p|0) = (0|p|0) + + (0|p|0)_, 



where 



(0|p|0) 



1 



+ 



2vrZ 



(0|p|0)_ ee 



2tt/ 



n=0 
j oo 



n=0 



2 E^ + ^-^) 

n=0 Z 



(25) 

(26) 
(27) 
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"i r 



- 1 r 

W c -[l 



"I 1 — 



-> b 



-2- 
-3- 
-4- 

-5- 



Figure 1: Graphical solution of eq.(30). The y axis is labeled in units of po = |^ 

We reconstruct V^// from (0|p|0) using 

K// = -(0|p|0)A -^6 2 , (28) 

where the second term is the tree level contribution. We note that in the absence of 
fluctuating electrmagnetic fields, the role of A is played by — in (0|p|0)+ and by 
in (0|p|0)_. Thus 

V eff (b) = 1ft ((0\p\0} + - (0|p|0)_) - ^6 2 . (29) 

We wish to look for extrema of the above away from the points b = ±[(2n + 1)^W — 2/i], 
since at these points the extrema picks up singular delta function contributions. The 
extrema satisfy 

(0|p|0) + -(0|p|0)_ = -6. (30) 

9 

Since the l.h.s. of eq.(^) is a complicated function of b, it is solved graphically (see 
fig.(l)). It is seen that if the external magnetic field is strong enough to satisfy B 2 > ^-§-, 
at least one non-zero solution exists which minimizes the effective potential. One such 
solution, which is for the situation where B 2 ~ ^-§-, is bo = u c . This, as has been 
discussed earlier in this section, leads to an effective Zeeman interaction term ^ij)* agij) . 
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Thus, we have the rather interesting result that a four fermi term in a strong external 
magnetic field can lead to an effective Zeeman interaction for the fermions. 

5. Relativistic Q.E.D. 

Having demonstrated the efficacy of our technique for the nonrelativistic case, we 
now go on to discuss the relativistic scenario. Although this is deemed unnecessary for 
most condensed matter oriented applications, there are situations where the fermionic 
dispersion relation may be taken to be approximately linear (e.g. in the neighbourhood 
of the Fermi surface) and our subsequent calculations would find applicability in pre- 
cisely these cases. Furthermore, as has been discussed at length in the introduction, the 
astrophysical and cosmological ramifications of this problem, namely fermions interact- 
ing through a gauge field, Abelian or otherwise, in a strong magnetic background, are 
manifold. The strong magnetic field can induce fundamental changes not only in the 
vacuum structure of such theories but also in the other sectors of the Fock space. 

Admittedly, in this work, we have addressed the problem of only planar fermions with 
an Abelian interaction. However, adding on the third spatial dimension entails only very 
minor modifications, namely, the single particle energies also come to depend quadrati- 
cally on the z-component of the momentum. This can be incorporated quite easily within 
our computational framework as can the extension to non-Abelian interactions, which 
actually bring in little more than notational complexity. 

Having enumerated these motivations, we now proceed to the actual computation. 
The generating functional is then given by 



where [i is the chemical potential and we use an euclidean metric with signature —1. The 
fermionic degrees of freedom can then be integrated out to yield the formal expression 
of the effective action as 




(31) 




Trlog^-m-^7 ] 



(32) 



Define a = a — — , and a 



a. Then we can express the effective action as 




Trlog[$>-m] , 



(33) 



with 



D 



(34) 



12 



Note that the commutator [ni,n 2 ] = -p. 

Upon functionally differentiating with respect to the gauge fields, the corresponding 
current is obtained as 

W= 2 Tr(x|{ ^'^ -}\*) ( 35 ) 
The derivative expansion can now be performed in an identical manner to the nonrela- 
tivistic case. The bra and ket are translated to the spacetime origin, and the operators 
to x^ + x M . We then formally expand all functions of x^ + x^ around x^. Also, the 
unitary transformations U, V, and W on the current operator are identical to those in 
the non-relativistic case. Upon performing these transformations the current takes the 
form: 



^ = >^>"W^-V m (36) 



with 



and 



0° = iexid'a + -XiXjd'&a + iex Xi& of , a° + ... (37) 



O l = ^Xi&a 1 + ^XiXjd^a 1 + iex d°a l + iexoXidP&a 1 + ^ x 2 (d°) 2 a l (38) 

Z o _ 

Again, we can use perturbation theory to expand : 

1 1 1^1 



+ ... (39) 



ifl— m — (p ifl— m itfl— m ifl— m 
Also, the algebraic relations between 7 matrices (see appendix A) can be used to express: 

1 1 



ifi- m IT 2 - m 2 + ^7 

Defining 



T-^+m) (40) 



S = U 2 -m 2 (41) 
M = (U 2 -m 2 ) 2 -^ 



Qi 
Q2 



_S_ 

Mo 
1 
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we obtain 

= (Qi + 7oQ 2 )(iJ&+ m) (42) 

and therefore 

(J M (x)) = -Tr(0\{ lfl ,(Q 1+lo Q 2 )(tfl+m) (43) 
+ (Qi + 7oQ 2 )(^+ m)0(Qi + 7oQ 2 )(^+ m)}|0> + . . . 



In order to calculate expectation values in the state |0), it will be necessary to ex- 
pand this state in the eigenbasis of the operator II 2 . It is instructive to compute the 
eigenspectrum of II 2 explicitly since this calculation will also motivate the use of guiding 
center coordinates . 

Consider the eigenvalue equation 

(II 2 - m 2 )^ = A^ (44) 

i.e. 

[(po + ¥f +P 2 i + (P2 - eB Xl ) 2 ]^ = -(A + m 2 )^ (45) 
One can solve this by separation of variables: 

x 2 , x ) = e~ lulx0 e~ ieBXx ' 2 <5>( Xl ) (46) 

This leads to the single differential equation for <3>: 

[-d\ + (eBfix 1 - X) 2 ]$ = -[A + m 2 + (u - ^) 2 ]$ (47) 

This is just the Schrodinger equation for a harmonic oscillator in one dimension., with 
eigenvalue spectrum (2n + l)// 2 . The full eigenvalue spectrum of the operator II 2 — m 2 
is therefore 

A n = -^±^-m 2 -(c,-^) 2 (48) 



The corresponding eigenstates may be labelled by the quantum numbers n, uj and X. It 
follows that 

SIX ,n,cu) — X n \X, n, to) (49) 
M \X,n,u) = r n \X,n,u) (50) 
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where r n = A 2 — h Therefore, 



Qx\X,n,u) = =P-\X,n,u) (51) 



r 



Q 2 \X,n,uj) = -^—\X,n,u) 



n 



and the completeness relation takes the form 

V /— dX\X,n,u){X,n,u\ = 1 (52) 
n J 2n 

We now define the operators 

X = -xx- l 2 U 2 = -l 2 p 2 (53) 

Y = -x 2 + fUx = -x 2 + l 2 v\ 

They have the properties 

X\X,n,u) = X\X,n,u) (54) 
[X,Y] = il 2 
X|0) = -z 2 n 2 |o> 
y\o) = / 2 nx|o) 

As we shall later see, these properties give us a powerful method of computing expectation 
values of quantities in the state |0). 

Consider eq. (|44|) . Since the trace operation is independent of the ordering of operators 
in the commutator, the expression for the current to first order in O is given by: 

(J M > = *eTr(0| 7At (Qi + ToQ 2 )(^+ m)0(Qi + ToQ 2 )(^+ m)|0) (55) 

This gives rise to 16 separate terms. The calculation is done in the following manner: 
First, in each term, the trace is taken with respect to the 7 matrices. This leaves products 
of Q's ,11's and O in each term. The dependence on the x^ in O is removed through the 
use of properties eq. fl53|) of the guiding center coordinates, as will be shown in a specific 
example below. Then the state |0) is expanded in the eigenstates of Q\. Remembering 
that II12 act as combinations of raising and lowering operators on these eigenstates, the 
expectation value can be computed in a straightforward, albeit tedious, manner. 
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To illustrate, consider the term — ie(0\ r YoQif]O r yoQ2^0) in (J ). Upon taking the 
trace with respect to the 7 matrices, this splits into 6 separate terms: 



2e(0|Q 1 n 1 O°Q 2 n 2 |0) - 2e(0|Q 1 n 2 O°Q 2 n 1 |0) - 2e(0|Q 1 n O 1 Q 2 n 2 |0) 
+2e(0|Q 1 n O 2 Q 2 n 1 |0) + 2e(0|Q 1 n 1 O 2 Q 2 n |0) - 2e(0|Q 1 n 2 O 1 Q 2 n |0) 



(56) 



We would like to obtain a contribution proportional to f 12 from these terms as a demon- 
stration. We note that contributions to f 12 are obtained from 0\ and 2 : 



ze 



0l - ? x 2/ - 



Oo 



/e 
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Thus the terms from eq.(|56"D which contribute to f 12 are: 



ze 2 [(o|Q 1 x 2 g 2 n 2 n |o)+(o|Q 1 n 2 x 2 Q 2 n |o)+(o|g 1 x 1 g 2 n 1 n |o^ 

Using eq.(|54D and the commutators in appendix B we get 



[r ~ ,s \s,ni}^nifiWo>-<o|^ 



12 



Now there is another term in (Jo) with Q\ — > Q 2 . This contributes to f 12 in the form: 



S 



S 



-^Kol^nj-^n^olo) + (oi-^— {s,n x }iio|o) - (o|— n 2 — n |o>]/ 



12 



l 2 iri M L "'" iJ ¥ 2 
Consider the simplest term: 



M 2 'Mo 



M L M< 







4e2 
I 5 " 

l 2 " 



/ 12 (o|n 1T - ^riolo) 



Mo 2 



(57) 



(0|n, X, w) (n|ni|ni) HriiK) (n 2 , X, w|0) 



n,ni,n2 



r 2 

ni 



The integration over X now gives po so we get 



2e 2 , 12 ^ r efoo, . , (n|IIi|rii)(ni|ni|n) 



n\,n 



2vr 



r 2 

ni 



(5f 



n + 1 n 

T 75 



r 2 

n+l 1 n—1 
1 



1 

_ L n+l i n. 
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simplifying the other terms in the same way we get 
2e2 xl2 f dw , \ 



2e% 12 ^ fdw f _ , A (n +l)(A„ + A n+1 ) (\ n | A n+1 \ 

^o- 7 r n r n+1 yr n r n+1 y 

The coefficient of the term proportional to f 12 in (jo) also receives contributions propor- 
tional to m. These are calculated to be 

2ie 2 



/12 g /• M (w + 1) ( A » + A "+i) ^ Mjn+i - g) + w A » - _ (2w + (59) 

J 2?r r 7 r 



A further simplification is in order in the above expressions. Let us define 

d n = — ^ -m 2 - (w - ifi) 2 . 

Thus 

1 1 / 2 / 1 1 \ 



(60) 



using this technique, any denominator containing arbitrary products of the Ts can be 
simplified to powers of a given d. Further, the numerators can trivially rewritten in terms 
of d since A n = d n — p. 

The final expression for (J ) and (Ji) can be expressed entirely in terms of w , and 
d n . It is quite straightforward to show, as in the nonrelativistic case, that the terms that 
potentially violate gauge invariance vanish identically. Compiling all the terms together 
we get: 

(J ) = C f 12 + C 1 d°d 1 a 1 + C 2 (d 1 ) 2 a° + ... (61) 
(Ji) = W + ^W + D 2 (9 C )V + ... 

where C , Cj and -D , -D« are constants to be specified later. Correspondingly, 

(n 00 (x,y)) = C 2 (9 1 ) 2 5(x-y)+C 2 (9 2 ) 2 5( a ;-y) + ... (62) 

(U Q1 (x,y)) = -C d 2 5(x-y) + C 1 d°d 1 5(x-y) + ... 

(U u (x,y)) = -C Q (d 2 ) 2 8(x-y)+D 2 (d°) 2 5(x-y) + ... 

(U 12 (x,y)) = -D d 2 d 1 8(x-y) + ... 

(n 02 (x,y)) = C d 1 S(x-y) + C 1 d°d 2 S(x-y) + ... 

(U 22 (x,y)) = -D (d 1 ) 2 5(x-y) + D 2 (d ) 2 5(x-y) + ... 
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Gauge invariance forces G\ = —C2 = D\ = — D2 = C so that the effective action that 
one gets after collecting all the terms and functionally integrating (11^) is 

S eff = ^J d 3 x£ 2 + ^-Jd 3 xB 2 + ^Jd 3 xa°B + ... (63) 

where 

p 2 OO p r-, 1 

D = D = — + 7n-- Q(y/m? + 2neB-\»\), (64) 



127T n=Q 



\/ m 2 + 2neB 



2w c 2. 

Ml 



and 



= ~ JHL-jjKL e(Vm 2 + 2ne£ - H) (65) 

32nw c ^ Vm 2 + 2neB 1 1 ' 



e 2 



g 71(10 ~ 1371) Q(vTO ^ 



H) 



167T ^ \/m 2 + 277ei? 

One sees that even in the relativistic case the coefficient of the euclidean S 2 vanishes for 
strong external magnetic field. The explicit dependence on the chemical potential \x can 
be removed by specifying the number of filled Landau levels. Physical quantities can be 
extracted from this action as in the nonrelativistic case. 

Discussions 

In this paper we have calculated the effective action of nonrelativistic and relativis- 
tic 2+1 dimensional QED in a strong external magnetic field using the inhomogeneity 
expansion method. We have explicitly addressed only the case of the parity odd model 
i.e. one fermion flavour , but it is straightforward to use our method in the parity even 



case by considering two fermion flavours with positive and negative mass fL5[ . For the 
nonrelativistic case we obtain a Chern-Simons like term in addition to the renormaliza- 
tion of the coefficients of S 2 and B 2 . In the relativistic case a true Chern-Simons term is 
generated and we calculate its coefficient as well as the leading term in the coefficients of 
E 2 and B 2 . For very stong external magnetic fields (I -C -\) we see that the coefficient of 
the euclidean S 2 term vanishes both in the relativistic and nonrelativistic case. This sug- 



gests that a dimensional reduction is indeed taking place [0. In this context, it would 
be interesting to study the effect of such a reduction in the case of Nambu Jona-Lasinio 
models in general and the 2+1 dimensional Thirring model in particular . 
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We have extracted just a few of the many possible physical quantities that can be 
obtained from these effective actions viz the dielectric constant and magnetic permeabil- 
ity. Considering a finite temperature variant we can further discuss the physics of 



plasmons in a magnetic field background. A particularly intriguing direction is to adapt 
these techniques to astrophysical and cosmological applications. 
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Appendix A. Dirac algebra and Trace relations 

In the relativistic calculation, the metric we use is g^ v = diag(— 1, —1,-1). One has 
to be careful, therefore, in raising and lowering indices. Every raising or lowering of an 
index picks up a negative sign. The canonical commutation relation takes the form 

[xn,Pu]=ig[iv (66) 

from which commutation relations involving EE's follow directly. 

The algebra of Dirac matrices in 2+1 dimensions, with a negative definite Euclidean 
metric, albeit straightforward, is not generally familiar. Here we quote the main relations 
which are used in the body of the paper. 

To start with, the Dirac matrices are defined as 71,2 = o"i,2, 7o = 0"3; where the cr's are 
the usual Pauli matrices, (note that 7* = — erj, because the metric has negative definite 
signature.) This leads to the relations 

IpTIv = -g^u + ie^plp (67) 

Hence 

{7m> 7^} = -2flV- ( 68 ) 

The trace relations are obtained as 

Tr7 M 7^ = -2^ (69) 
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Tr7 M 7^ la 7/375 = -^(g^apS + gafS^fiuS + g^S^va - gSa^^fs) 

viai pi P i a = -2(gnvgapgpa - gu a gvpg P a + g^gvagpa - gapgu P gva 

+gapgp^gyp — gp,vga P gi3a + g^vgaagpp — ^p,vp^afia + 



and so on. The trace of an odd number of 7 matrices does not vanish as it does in 
3+1 dimensions. 

Appendix B. Some useful relations 

First, for the non-relativistic case 



[n 2) — ] 

-no 
-no 

r. 1 , 



ml 2 H 'H 

i 



mPH 2 H 



H n 



1 

-"0 



For the relativistic case 

[xo, Qi\ 

[x , Q2] 

\fiiM 
[tli,Qi] 
[ri 2 ,Qi] 
[n 2 ,g 2 ] 



-4iS—Sfl — + 2ift — 
M °M M 

4 * 1 9ft 1 
"7^Mo Mio Mo 

~fraclM {S,il 2 }± 

2i S ~ 1 2i ~ 1 
TM~ |S ' n2| Mo + F n2 Mo 



2i 

2i 1 



1 1 1 

_L/5- ni}— -ni— 

Mo 1 J M M . 
1 



l^ {5,IIl} M 



(70) 



(71) 
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